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Abstract 

We combine magneto-optical imaging and a magnetic field pulse technique to study domain wall 
dynamics in a ferromagnetic (Ga,Mn)As layer with perpendicular easy axis. Contrary to ultrathin 
metallic layers, the depinning field is found to be smaller than the Walker field, thereby allowing for 
the observation of the steady and precessional flow regimes. The domain wall width and damping 
parameters are determined self-consistently. The damping, 30 times larger than the one deduced 
from ferromagnetic resonance, is shown to essentially originate from the non-conservation of the 
magnetization modulus. An unpredicted damping resonance and a dissipation regime associated 
with the existence of horizontal Bloch lines are also revealed. 

PACS numbers: 75.50.Pp, 75.60.Ch, 75.70.Ak 
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In ferromagnetic systems, domain wall (DW) motion driven by a magnetic field 



or a spin-polarized current 



flflfl 



, 13, |8| presents a variety of dynamical regimes. Depending 
on the field strength, several regimes characterized by the dynamics of the magnetization 
vector inside the DW and the DW mobility (field derivative of the velocity) are predicted to 
occur. Theoretically, the dissipation-limited regimes were mostly investigated in a system 
consisting of an ideal ferromagnetic film with uniaxial perpendicular anisotropy, subject to 
a magnetic field H parallel to the easy axis . In the Walker steady regime, the DW 

structure is stationary and a linear velocity v = \iH is expected up to the Walker field {2, 4|. 
Above this field the precession of the DW magnetization around H leads to a DW back and 
forth motion, reducing its average velocity |3(. This regime may become unstable, leading 
to the nucleation and propagation of Bloch lines inside the DW [l, In the high-field 
range of the precessional regime the damping torque becomes large enough to move the DW 
with linear velocity but reduced mobility with respect to the Walker regime. These various 



dynamical regimes could be observed in the past in micrometer-thick garnet films 
and very recently in nanowires with in-plane magnetization 1(3, However, in ultrathin 



metallic ferromagnetic films with perpendicular magnetization dissipation- limited regimes 
are masked by thermally-activated creep and depinning regimes jlj . . For instance, in 
Pt/Co/Pt layers only the high field linear precessional regime could be observed owing to 



the strong pinning of DWs 



12|. 



In ferromagnetic semiconductors, DW dynamics has been explored only recently. The 
creep regime has been observed in GaMnAs thin films and wires (si, [l^j]. However, the 
strength of the applied field was too small to reach the dissipation-limited regimes. The 
observation of these regimes is of prime importance to determine the relevant dissipation 
processes involved in DW motion. In particular, the nature of the ferromagnetism (hole- 
mediated interaction between diluted Mn ions Q, Q) could lead to specific features in the 
DW dynamics. Moreover, the understanding of dissipation processes in GaMnAs should 
have important implications for the study of current-driven DW motion. It should help 
to discriminate between the damping contribution and the spin transfer one, which is not 
yet well understood [7| . GaMnAs is a good candidate for experimental investigation of this 
question since the current density required to move a DW was found to be two orders of 
magnitude smaller than in metallic nanowires 

In this letter, we report on DW dynamics in GaMnAs over a wide range of magnetic 
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field and temperature values. A quantitative analysis of the DW dynamics leads to the 
identification of the steady and precessional regimes and to a self-consistent determination 
of the DW width and damping parameters. This damping is shown to involve the variation 
of the magnetization magnitude. We also identify a dynamical regime consistent with the 
existence of horizontal Bloch lines. Finally, DW dynamics reveals an unpredicted velocity 
peak. 

The sample consists of an annealed Gao.93Mno.07 As epilayer of thickness d=50 nm grown 
on a relaxed Gao.902 in o.o98As buffer deposited on a GaAs substrate. After annealing, the 
Curie temperature Tq is 130 K and the magnetic easy axis is perpendicular to the sample 
plane. Kerr microscopy is used for the direct observation of DW motion (see Ref. jisl . [19I for 
more details). The sample is placed in a helium-flow cryostat. The DW velocity is measured 



using a magnetic pulse field technique 
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20j. For low velocity (v < 10 3 m s 1 ), pulses are 



generated by a conventional external coil (rise time ~ 100 ms, maximum amplitude 60 mT). 
The DW velocity is determined from a set of snapshots tracking the DW position during 
a field plateau (see image (a) in Fig. [T]). For higher velocity (v > 10~ 3 m s _1 ), pulses are 
generated by a small coil (rise time ~ 200 ns, maximum amplitude 250 mT) of diameter 
« 1 mm placed inside the cryostat, onto the sample surface. DW motion is driven by a 
series of pulses of constant amplitude H and increasing duration r. The snapshots recorded 
before and after each pulse are used to determine the DW displacement as a function of r 
as shown in images (b) and (c) of Fig. [U The DW velocity reported in Fig. [T] corresponds 
to a linear fit of this curve (not shown). This procedure eliminates the effects of the pulse 
rise and decay times. 

A typical velocity curve for field-driven DW dynamics is shown in Fig. [1] (T = 80 K). 
Several dynamical regimes can be identified. For 1 mT< fioH < 2.5 mT - creep regime - DWs 
are rough and present spatially inhomogeneous displacements (image (a) in Fig. []]). The 
free DW motion is impeded by defects, some of them appearing as pointlike or long linear 
defects Q] ■ To determine the average velocity, DW displacements were only measured in 
areas with no strong pinning defects. The results are reported in the inset of Fig.[U v is found 
to increase over 6 orders of magnitude. The good agreement with a fit v = Voexp(—aH~ 1 / 4: ) 
suggests that DWs follow a creep regime described by the motion of an elastic string in the 



presence of a random pinning potential 



131 ] . However, the systematic investigation of this 



regime is out of the scope of this paper. For 2.5 mT< noH < 8 mT - depinning regime - the 
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velocity varies linearly with the field. The DW roughness decreases (image (b) in Fig. [T]) and 
the domains expand almost isotropically. One still notices pointlike defects around which the 
DWs fold up, creating lamellar domains Q- This regime, separating the creep and flow 
regimes, corresponds to the depinning regime. For fi H > 8 mT - flow regimes - the DWs 
are smooth and the displacements homogeneous (image (c) in Fig. [T]), thereby indicating 
that the DW dynamics is no more limited by pinning. The v(H) curve (Fig. [1]) presents the 
main features of the predicted dissipation-limited DW dynamics , |4| . A velocity peak 

(v = 10.5 m s _1 for [i H = 8.2 mT) is followed by a region with negative differential mobility 
(8.2 mT< fioH < 35 mT). The DW velocity decreases down to 7 m s _1 for n^H = 35 mT. 
For fi H > 50 mT, the velocity increases again with the field, v is proportional to H 
as expected for the high field precessional regime, except near fi H = 92 mT where an 
unpredicted second velocity peak is observed. In order to determine whether those features 
of the DW dynamics are systematically observed, measurements were performed over a wide 
temperature range (0.03< T/Tc <0.92). 

The results are reported in Fig. [2j Qualitatively, the v (H) curves show a weak dependence 
on temperature for the depinning and flow regimes. The main difference concerns the upper 
boundary H up of the investigated field range. At low temperature (T = 4-50 K), H up 
corresponds to the maximum available field amplitude (250 mT). At higher temperature 
(T = 65-120 K) H up is limited by nucleation. The distance between domains nucleated 
during a pulse is too small for an accurate measurement of the DW displacement. The 
anomalous velocity peak in the high field precessional regime is systematically observed, as 
shown by the arrows in Fig. [2j 

In the following the main features of the DW dynamics are analyzed quantitatively. 
In uniaxial ferromagnetic films, the anisotropy is characterized by the quality factor 
Q = 2K u /ji M 2 , where K u is the uniaxial anisotropy constant and M s the saturation magne- 
tization [4]. Q is found in the range 8.6-14, which denotes strong uniaxial anisotropy [2l|. In 
that case, the velocity in the Walker steady regime is given by v st (H) = fj, st fj, H = •yAfi H/a, 
where 7 is the gyromagnetic factor (1.76 10 11 Hz T" 1 ), A the DW width parameter, and a 
the damping parameter [4J]. In the high field range of the precessional regime the velocity is 
given by v prec (H) = /j prec fx H = 7A/i #a/ (1 + a 2 ) [4|. 

Those predictions are compared to the observed dynamical regimes. Since the flow regime 
is reached close to the velocity peak the mobility fi st is obtained by adjusting a straight line 
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tangent to the experimental curve (Fig. [3] (a) and (b)). /i pr ec is obtained from a linear fit of 
the high-field regime (Fig. E^a)). The ratio fi pre c/ 'fJ-st = a V(l + ct 2 ) yields a. A is obtained 
as A = fi prec (l + a 2 )/cry. As shown in Fig. HJa) A varies weakly with the temperature. 
Moreover, the A(T) values lie between the two boundaries deduced independently from 
domain theory for the same sample 21(. This very good quantitative agreement confirms 
that the procedure used for the determination of A and a is relevant. This also demonstrates 
that the precessional and steady regimes are indeed observed experimentally, the latter being 
reached only near the Walker velocity peak as shown in Fig. [3]^b). 

The damping coefficient a ~ 0.3 is also weakly dependent on the temperature (Fig. IH(b)). 
A similar value (a ~ 0.15) was obtained from the decay of the magnetization precession in- 



duced by optical excitation 



221 ] . Surprisingly, the DW damping is more than one order of 



magnitude larger than the damping deduced from the frequency dependence of the ferromag- 
netic resonance (FMR) linewidth: a ~ 0.01 for this sample 23], a value in agreement with 
theoretical predictions for GaMnAs (a ~ 0.02 — 0.03) 24j. A similar discrepancy was re- 



ported for garnet films with small a- values 



Theoretically, the DW equation of motion 



is derived from the Landau-Lifshitz-Gilbert (LLG) equation for the magnetization. The dis- 
sipation is classically described by the phenomenological Gilbert damping coefficient a that 
accounts for the relaxation of the direction of the magnetization vector , |4j. However, 

dissipation processes are expected to differ for uniform magnetization (FMR) and for a mov- 
ing DW. A deviation from the linear steady regime is predicted when the interaction of the 
moving DW with thermal magnons is taken into account [26]. However, this interaction can 
be safely discarded here since it would give a significant contribution only at high velocity, 
above ~ 300 m s _1 for this sample. In the low velocity range, additional relaxation terms in 
the LLG equation, taking into account the non-conservation of the magnetization modulus 
and heat exchange with a thermostat, can lead to a DW dynamical damping larger than the 
Gilbert (FMR) damping [27|, [28]. For a magnetization response time tm = X\\llM s otFMR 
shorter than the DW transit time r t = A/v, the dissipation for DW motion in the steady 



regime is described by a dynamical damping [27| aow = oifmr 



1 + 16 (4ttx||Q) /3a 



FMR 



where \\\ is the longitudinal magnetic susceptibility (CGS units). Using this equation with 
a fmr ~ 0.01 yields \\\ ~ 10~ 4 , in reasonable agreement with theoretical estimations based 
on spin-wave theory j3; 13 • Taking into account the anisotropy gap, one finds x\\ m the 
range 4 10~ 6 -10~ 5 for T = 12 K and 8 10~ 5 -2 10" 4 for T = 80 K. For the measured velocities 
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one finds tm <C r t , thereby justifying the use of this model. We point out that the theoretical 
prediction of enhanced DW dynamical damping in the steady regime also accounts well for 
the damping in the high- field precessional regime. 

Let us now extend the analysis of the experimental velocity curves beyond the ID theory 
of DW motion. As calculated and shown in Fig. E](b), the ID theory predicts a maximum 
velocity (Walker velocity vw = lAfi Q M s /2) at the Walker field Hw = otM s /2. For H just 
above Hw the DW back and forth motion due to the DW magnetization precession around 
the applied field should lead to a decrease of the time-averaged velocity and hence to a 
region of negative differential mobility 3]. In contrast, the experimental v (H) curve just 
shows a change of slope with still a constant positive mobility, yet strongly reduced with 
respect to the steady regime (Fig. [3(b)). In this field range, the DW structure is expected to 
be unstable. A solution for DW propagation with generation and propagation of horizontal 
Bloch lines through the film has been proposed For sufficiently large a this model 
predicts a viscous-like drag with decreased mobility hbl = 7 A/ [a (1 + n 2 A/2a 2 a)} with 
A = Ay/Q the exchange length and a of the order of the film thickness The linear fit of 
the velocity curve just above the Walker field for the set of investigated temperatures yields 
a in the range 86±50 nm, consistent with the sample thickness d — 50 nm. Given the fact 
that the exchange length d/10 < A < d/5 is not very much smaller than d, contrary to the 
assumption of the model, the agreement is quite satisfactory. It strongly suggests that DW 
motion above the Walker field is slowed down by the repetitive generation of one Bloch line 
at one surface of the film and subsequent propagation and annihilation at the other surface. 

Let us now discuss the intriguing velocity peak observed around 90-120 mT (Fig. [2]). To 
our knowledge, such a peak has been neither predicted nor observed in ferromagnetic sys- 
tems. Since it occurs within the linear precessional regime, where the velocity is proportional 
to the damping parameter, it can be ascribed to a damping resonance. In order to char- 
acterize its temperature dependence, the velocity curves are fitted using a field-dependent 
damping a(H) = a + 5aexp — (H — H p ) 2 /2a 2 . As shown in Fig. [2] (inset) the resonance 
amplitude (<5a=0.04 at 60 K) decreases linearly with the temperature, extrapolating to zero 
at T 120 K, close to Tq- The resonance field H p varies slightly with the temperature. The 



corresponding energy Too = Ti^Hq^J H 2 — H w /(l+a 2 ) is of the order of 10 fieV (u = 17 GHz). 
It may correspond to transitions between confined states of volume magnons, as calculated 
from the dispersion curve using the spin stiffness constant determined from this work and 
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ref. 



21] . The damping resonance might also be related to DW excitations (flexural modes), 



whose energy spectrum lies inside the anisotropy gap of the volume magnons (43 /xeV at 
80 K). These excitations have been calculated for a static or moving DW in the steady 
regime j^, Q] but not in the precessional one. 

Despite the small value of the saturation magnetization and hence of the Walker field 
in GaMnAs, we could observe the dissipation-limited flow regimes beyond the creep and 
depinning regimes. The steady as well as the precessional regime are described by a DW 
dynamical damping, which is found to be much larger than the FMR damping, in agreement 
with theoretical predictions considering non-conservation of the magnetization modulus. 
However, contrary to the classical assumption, a single, field-independent, damping constant 
cannot account for the whole DW dynamics. The existence of a damping resonance inside 
the precessional regime suggests an additional dissipation mechanism and calls for further 
theoretical investigations of dissipation processes in this regime. 
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FIGURE CAPTIONS 



FIG. 1 Magnetic field dependence of the DW velocity at T = 80 K. The error bars 



d region in 



131 ] . Labels 



represent the standard deviation of the velocity distribution. Inset: low fie _ 
semi-logarithmic scale with a fit (dashed line) according to the creep model 
(a), (b), and (c) refer to the corresponding images in the creep, depinning and flow regimes, 
respectively. Images (b) and (c) are differential images showing the DW displacement as 
indicated by black arrows. 



FIG. 2 (Color online) Temperature dependence of the DW field- velocity curves. Each 
curve is up-shifted by 10 m s -1 with respect to the previous one. Inset: temperature de- 
pendence of the resonance field H p and damping da of the velocity peak indicated by arrows. 



FIG. 3 (Color online) Comparison between theoretical predictions and experimental 
results (black circles). Dashed blue lines in (a): linear velocity in the steady and precessional 
regimes with mobility fi st and n preC) respectively. Solid blue curves in (b): velocity in the 
steady regime (below H w ) and in the unstable precessional regime (above Hw) calculated 
using the experimentally determined values of A and a (this work) and M s [2lj; dotted red 
line: linear fit of the velocity in the Bloch line regime. 



FIG. 4 (Color online) (a) Comparison between the DW parameter A obtained from 
DW dynamics (full circles) and the boundaries for A obtained from domain theory (open 
symbols) 21]. (b) Temperature dependence of the DW damping parameter a. 
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